Quantum algorithms which accept hot qubit inputs 
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Realistic physical implementations of quantum computers can entail tradeoffs which depart from 
the ideal model of quantum computation. Although these tradeoffs have allowed successful demon- 
stration of certain quantum algorithms, a crucial question is whether they fundamentally limit the 
conrputational capacity of such machines. We study the limitations of a quantunr computation 
model in which only ensemble averages of measurement observables are accessible. Furthermore, 
we stipulate that input qubits may only be prepared in highly random, "hot" mixed states. In 
general, these limitations are believed to dramatically detract from the computational power of the 
system. However, we construct a class of algorithms for this limited model, which, surprisingly, 
are polynomially equivalent to the ideal case. This class includes the well known Deutsch-Jozsa 
algorithm. 



theless, it is still possible that new techniques will be 
found which make them equivalent; this has not been 
proven impossible, but such a result would be unex- 
pected. Second, the limitations of the BQC model are 
not fundamental - they arise from practical considera- 
tions which significantly simplify experimental realiza- 
tion. Might there be other kinds of restrictions to the 
SQC model which apply to the most important quan- 
tum algorithms, and also simplify implementation while 
providing a provably polynomially equivalent model of 
quantum computation? 

In this Letter, we prove the existence of non-trivial 
overlap between the BQC (with thermal inputs) and SQC 
models: a surprising subclass of quantum algorithms for 
which BQC and SQC are polynomially equivalent. This 
includes the well known Deutsch-Jozsa (DJ) algorithm. 
The equivalence arises from the robustness of the algo- 
rithms to independent bit flip errors which occur before 
the computation. Our result does not resolve the general 
question of equivalence between BQC and SQC, but it 
provides new insight into the degrees of freedom which 
quantum algorithms allow different computational mod- 
els. 

It is convenient to begin with the Deutsch-Jozsa algo- 
ithm M (as improved in Jig] ), which solves the follow- 
ing problem. Alice has a function / : {0, 1}™ — ► {0, 1} 
which is unknown to Bob. She promises that / is ei- 
ther constant or balanced (f(x) = for exactly half of 
the inputs and f(x) = 1 otherwise). Bob wants to find 
out whether / is constant or balanced with the minimum 
number of queries. With a deterministic classical com- 
puter, 2™ _1 + 1 queries are required in the worst case 
to determine the type of the function. However, in the 
SQC model, the following algorithm can be performed: 
1. Starting from the state |0)®"|1), Bob prepares a query 
input register in a superposition of all basis states using 
the Hadamard transformation H, and an ancilla qubit 



The discovery of fast quantum algorithms [pJ-W with 
no classical counterparts has led to an increasing interest 
in the experimental realization of quantum computers. 
However, in practice, it is expected that realistic imple- 
mentations of these machines will suffer from imperfec- 
tions which depart from the ideal model of quantum com- 
putation. Ideally, standard quantum computation starts 
with a system prepared in a pure state such as |00 • • • 0), 
and unitary transforms are applied systematically as pre- 
scribed by the algorithm, obtaining some final state \ip). 
When measured, this state collapses with high proba- 
bility into a small number of possible outcomes, which 
reveal the desired answer. 

In contrast, recent demonstrations of simple quantum 
algorithms M pf have employed a far different model of 
quantum computation, in which the system is a bulk en- 
semble of a large (but finite) number of quantum com- 
puters operating in parallel, prepared initially in a highly 
mixed state, and with the only accessible measurement 
result being an average of the observable over the en- 
semble. Despite these stringent limitations, it has been 
shown how quantum computation can indeed be per- 
formed, by creating "effective pure states" [[io|-|l3|], and 
by modifying quantum algorithms to produce determin- 
istic results which do not average away p0[ . However, it 
is understood that using present techniques, this ability 
comes at a cost: either an exponential reduction in signal 
strength as a function of the size of the computer, or a 
linear reduction in the number of usable qubits E3] , with 
all other resources held constant. 

These observations lead to two important points mo- 
tivating the present work. First, it is an open question 
whether or not the limited "bulk quantum computer" 
(BQC) model described above is as powerful as the stan- 
dard quantum computer (SQC) model. Previous work 
has shown that BQC with thermal inputs is strictly less 
powerful than SQC in the presence of oracles [15[. Non- 



in the state H\l) = (|0) - |l))/\/2. Bob sends both the 
input and ancilla to Alice. 2. Alice evaluates f(x), stores 
the answer in the ancilla qubit and sends the qubits back 
to Bob afterwards. 3. Bob applies H to the input reg- 
ister and the ancilla qubit. These three steps can be 
represented by the following expressions: 
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where all summations are over {0, 1}™. The final state of 
the input register is therefore 
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If / is constant, \g(y)\ 2 = 5(y,0); if / is balanced, 
|g(0)| 2 = 0. When Bob projects the qubits along the 
basis states, he obtains a definite value of y to decide 
whether / is constant or balanced. The problem is there- 
fore solved with only one query in the SQC model. 

Does the DJ algorithm also work efficiently in the BQC 
model? First, let us consider a less restrictive model than 
BQC, which we shall term BQCp. This model is the 
same as BQC, but with the initial state of the quantum 
computers prepared in the pure state 1 00 • • • 0). For this 
model, a resolution to our question comes from looking 
at the DJ algorithm from the following, different point of 
view. 

Projection onto the basis states, which are eigenstates 
oia z (a z \x) = (2x — l)\x),x G {0, 1}), performs mesaure- 
ments of the eigenvalues ±1 from each qubit. These 
eigenvalues are —1 for all qubits when / is constant and 
+1 for at least one qubit when / is balanced. Thus the 
DJ algorithm can be viewed as distingushing these two 
possible outcomes. For BQCp, the output from the i th 
qubit is the expectation value of o~ Zi , given by 
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where yi is the i th bit of y. If f(x) is constant, £^ = 

— 1 for all i. If / is balanced, 3i such that E i ® p > — 1. 
Due to the ensemble average over measurement results 
required in BQCp, the two types of functions may not 



be distinguishable, as they are in SQC. The principle 
question is, can the two cases be distinguished without 
exponentially increasing the space-time complexity? 

We answer in the affirmative by analyzing the following 
measure of distinguishibility, defined as 
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That is, e(n) is the difference in the signals from the con- 
stant function and the worst balanced function Alice can 
apply. We show that BQCp and SQC are polynomially 
equivalent by proving the following fact that e(n) is lower 
bounded by the inverse of a polynomial in n. 
Fact 1: Vra, e(n) > 2/n. 

Proof. Since E i ^ on g t = — 1 is independent of / and the 
maximum over i is lower bounded by the average, from 
Eq.(0), 
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Using Eq. 
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where the crucial fact |<?(0)| 2 = inherent in DJ is used 
to omit the y = term in Eq. ([LOh. We have also used 



the normalization J2 v =o \d(y)\ 2 — 1- Since Vy > 1, 
Yl7=i Vi — ■"•) ^ follows that 
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The desired result is obtained by substituting Eq.(|l] 
into Eq.(|). D 
Since e(n) > 2/n, it can be classically amplified to be- 
come detectable in time 0(n 2 ). Thus, BQCp can solve 
the DJ problem with 0(n 2 ) queries, exponentially better 
than deterministic classical computation. 

Now we return to the BQC model, which is more re- 
alistic than BQCp, as thermal initial states are typically 
more experimentally feasible to prepare than pure states. 
The thermal initial state of BQC can be conveniently ex- 
pressed in the notation of |L8| as [in) = ®ky/Pk\k), where 
[•) denotes a probabilistic mixture over some ensemble 
of pure states, \k) = \ki,k 2 , . . . , k n ) with ki £ {0,1}, 
and pk is the probability of having the pure state |fc). 
When the i th thermal qubit being |0) has probability qi, 
Vk = ^hli~ '(1 — 1i) ki represents an uncorrelated initial 
distribution. Applying the DJ algorithm to the thermal 
state directly, one obtains 
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where Uf denotes the restriction of Uf to the input reg- 
ister. For simplicity, we have assumed that the ancilla 
qubit is pure. It will be clear why it does not affect our 
conclusion. We prove in the following that BQC and 
BQCp are polynomially equivalent for the DJ algorithm. 
Fact 2: For the DJ algorithm, all BQC measurement re- 
sults are proportional to those from BQCp with constant 
proportionalities independent of n. 
Proof: First, we note that for each qubit, \ki) = cr^|0), 



therefore \k) = \k\, . 
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• 0) , which 
is abbreviated as cr x \0). Second, from Eq. (||), it is clear 
that Uf only changes the relative signs of the states \x); 
hence, Uf commutes with o~ z on any input qubit. It fol- 
lows that [a£, H® n U f H® n ] = H^ n [a^, U f ]H® n = 0. The 
output state of the DJ algorithm in BQC as given by 
Eq. (O) can therefore be written as: 
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Hence, the signal from the i qubit is given by 
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where we have used efe — Pr(fc^ = 0). □ 
From Facts 1 and 2, it follows that e(n) is at least 
(2q' — 1)- for q' = min(qi, ■ ■ ■ ,q n )- The fact that the 
ancilla is a thermal state will only reduce e(n) by a con- 
stant factor of two. This establishes that the DJ algo- 
rithm operating with the BQC model also requires only 
0(n 2 ) queries, produces an output with no loss in signal 
strength, and achieves an exponential speed up relative 
to deterministic classical computation. 

This surprising result can be generalized to define 
an entire class of quantum algorithms which are poly- 
nomially equivalent on both the BQC and SQC mod- 
els. We shall refer to these as "hot qubit algorithms" 
(HQA). Let the initial state of the BQC be n uncorre- 
lated thermal qubits. This state is related to p B Q Cp = 
|00 • • • 0) (00 • • • 0| through an independent bit flip er- 
ror operation £ such that p B Q c = £(p Bt ^ Cp ) 1 where 
£{p) = (aJU&Xp) = EkPkVxPVx with £i(p) = q lP + 
(1 — qi)o~ Xi po~ Xi . Let the measurement signal from the 



i th qubit be Ei{p) = (cr Zi ) = tr(pa Zi ). A quantum algo- 
rithm U is an HQA if it satisfies the following two condi- 
tions. (1) BQCp can implement U with at most polyno- 
mial slowdown compared with using the SQC model. (2) 
Ei{J-{p)) = CiEi(p) where Cj are independent of n and T 
is defined byUo£ — Toll with U{p) = UpU^. Note, when 
U and £ are given, T is determined by J-(p) = J2 k B^pB\ 
with B k = y/pffJo*Tfl. 

We show here that condition (2) ensures BQC and 
BQCp are polynomially equivalent. Since the initial 
state of the BQC is p BQC = £(p BQCp ), implemenent- 
ing IA in BQC is equivalent to implementing hi o £ in 
BQCp. By condition (2), U{p BQG ) = ToU{p BQCp ) and 
Ei\u\p BQC )) = c i E l (U(p BQCp ). In other words, under 
condition (2), the initial state imperfection only causes 
a constant signal loss without changing the output infor- 
mation, hence, polynomially relating BQC and BQCp. 
Together with condition (1), BQC and SQC are polyno- 
mially related in implementing HQA. 

We now consider what algorithms are HQAs using the 
criteria we have developed. From previous discussions, 
it follows that the DJ algorithm is an HQA with T = £ 
and Ei(T(p)) = (2^ — l)Ei(p). Another problem solved 
by an HQA is the parity problem |16| , which determines 
an unknown y G {0, 1}™ using the minimum number of 
queries of the inner product function / : {0, 1}™ — > {0, 1}, 
f(x) = x ■ y. This problem is solvable by the same pro- 
cedure given by Eqs. (0)-(Q). Now with f(x) = x ■ y, the 
final state following from Eq. (J3|) is, 
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BQC outputs Ef QC = (2 ft - l)E? QCp = (2 ft - l) Vl . 
Note that in contrast to the DJ problem, the parity 
problem has a deterministic answer such that no clas- 
sical amplification is needed. Since the reduction in the 
signal, 2(ft— 1, is independent of n, there is sufficient confi- 
dence to determine y with only one query, the same as for 
the SQC model. Classical information theory, however, 
shows at least n queries are needed when y is randomly 
distributed. 

Condition (1) for an algorithm to be an HQA may seem 
to be overly restrictive, but in reality it is not. We have 
found efficient modifications of all known quantum algo- 
rithms |16|,[l7j f° r the BQCp model. The general method 
for doing this, as pointed out originally in JTO] , is to per- 
form the usual classical post-computation (such as con- 
tinued fraction expansions, and inner products) as part 
of the quantum algorithm, on each quantum computer 
in the ensemble. Furthermore, since these quantum al- 
gorithms generally produce solutions which are efficiently 
verifiable, this procedure in effect determinizes the quan- 
tum algorithm and translates probabilistic operation into 
a bounded constant signal loss. 



Nevertheless, the importance of condition (1) for an 
algorithm to be an HQA is paramount, due to the dif- 
ficulty to satisfy both conditions (1) and (2) simultane- 
ously. For example, Simon's algorithm H, unmodified, 
satisfies condition (2). However, BQCp outputs a cor- 
rect answer only for a few cases, but otherwise erases 
all useful information on the answer. Therefore, BQCp 
fails to solve Simon's problem with the original algorithm 
even though BQCp and BQC are polynomially related. 
Meanwhile, the modified Simon's algorithm for BQCp 
does not satisfy condition (2). Thus, we have not yet 
been able to construct an HQA for Simon's problem. 

HQA is easily generalized to computation models in 
which the elementary carriers of quantum information 
are q-sny states for q being a prime power. In this ex- 
tension, the bit flip o x is replaced by the generalized 



Pauli operator X q defined as 



X q \i) = |i + l), where 
Then, £ relat- 
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addition is taken to be modulo q. 
ing the BQC initial state p BQC to |00 • ■ • 0) can sim- 
ilarly be constructed. When implementing the algo- 
rithm U, BQC and BQCp are polynomially related pro- 
vided the corresponding operation T satisfies E i (J r (p)) = 
CiEi(p), with Ci independent of n. An algorithm satisfy- 
ing this is U 

denotes the discrete Fourier Transform, 
~T S«=o e l27rx ' y / q \y), and Uf commutes with the phase 
operator Z q , generalized similar to X q . 

In conclusion, we have shown the existence of a class 
of algorithms, HQA, which are polynomially equivalently 
realizable with the standard quantum computer model, 
and the bulk quantum computer model with either pure 
or thermally mixed inputs. These algorithms can thus re- 
tain many advantages provided by quantum algorithms 
over the classical, even when implemented on physical 
systems with limited capabilities. In particular, for non- 
trivial algorithms in HQA such as the DJ algorithm, 
which require entangling unitary operations (since sin- 
gle qubit operations can only implement 2(2" — 1) bal- 
anced functions out of ( 2 «-i) m total), the necessity of 
such nonclassical operations implies that classical simu- 
lation of BQC would be difficult even though BQC might 
involves only separable states [|19|. HQA can be gener- 
alized in several ways. First, the BQC initial state can 
be any state, related to |00 - - ■ 0) through a general er- 
ror operation £, other than the thermal state. Second, 
modifications to the original algorithm and non-unitary 
initialization can be performed as long as they do not 
detract from the complexity of the original algorithm U. 
More concretely, let U satisfy condition (1). If there ex- 
ists an operation V, possibly non-unitary, such that !F, 
defined by Vo£ = ToU, gives E i (J 7 (p)) = CiEi(p) with a 
independent of n, then U is an HQA. The only constraint 
is that the implementations of V and ti have polynomi- 
ally equivalent complexity. These results indicate that 



a deep relationship exists between hot qubit algorithms 
and algorithms which do not spread initial errors, or al- 
gorithms involving quantum error correction. It is likely 
that a good place to look for applications of the BQC 
model would be in quantum algorithms which are robust 
against the initial uncorrelated bit flip errors. Further 
study of such relationships may lead to other interest- 
ing quantum algorithms for bulk quantum computation 
with imperfect input states, and other physical models of 
quantum computation which are realizable but limited by 
practical considerations. 
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